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ABSTRACT 

Context. Superluminous type la supernovae (SNe la) may be explained by super-Chandrasekhar-mass explosions of rapidly rotating white 
dwarfs (WDs). In a preceding paper, we showed that the deflagration scenario applied to rapidly rotating WDs generates explosions that cannot 
explain the majority of SNe la. 

Aims. Rotation of the progenitor star allows super-Chandrasekhar mass WDs to form that have a shallower density stratification. We use simple 
estimates of the production of intermediate and iron group elements in pure detonations of rapidly rotating WDs to assess their viability in 
explaining rare SNe la. 

Methods. We numerically construct WDs in hydrostatic equilibrium that rotate according to a variety of rotation laws. The explosion products 
are estimated by considering the density stratification and by evaluating the result of hydrodynamics simulations. 

Results. We show that a significant amount of intermediate mass elements is produced for theoretically motivated rotation laws, even for 
prompt detonations of WDs. 

Conclusions. Rapidly rotating WDs that detonate may provide an explanation of rare superluminous SNe la in terms of both burning species 
and explosion kinematics. 

Key words. Stars: supernovae: general - Hydrodynamics - Methods: numerical 



1. Introduction 

Type la supernovae (SNe la) were long believed to form a rela- 
tively homogeneous class of events in terms of their spectra, 
peak luminosities, and light curves. Surveys, however, have 
shown that the distribution of SNe la properties is substantially 
broader than previously anticipated, including those of h ighly 
peculiar events such as SN 20 07ax dKasliwal et alj|2008l). SN 



ics as well as spectra, light curves, and nucleosynthesis of 
type la supernovae. So far, a major disadvantage of this 
model has been that the conditions for a deflagration-to- 
deto nation transition (DDT) as proposed by iKhokhlovl (I1991I) 
and IWoosley & Weaven d 19941). are qu estionable in the con 



text o f WD matter dNiemeyen [1999). Although IPan et al 



2005hk dStanishevetalJ 120071) . and SN 2003fg dHowell et al 
2006). The latter, in particular, has been interpreted in terms 



of a super-Chandra sekhar-mass explos ion of a rapidly rotating 
white dwarf (WD) djefferv et al.ll2006l) . 

Motivated by the goal of explaining the presumed ho- 
mogeneity of SNe la within a single explosion scenario, 
most multi-dimensional hydrodynamic al simulations so far 



(2008) proposed a theoretical ex planation of de l ayed deto- 
nations, the n umerical studies by I Wooslev et al.1 ([2008) and 
Schmidt et al. (2009) nevertheless impose severe constraints on 
DDTs. Prompt detonations, on the other hand, have been con- 
sidered to be in feasible for ex plaining SNe la since the pio- 
neering work of Arnett dl969h . as too little material is burned 
at sufficiently low densities to produce intermediate-mass ele- 
ments (IMEs). 



2005 



2006: lRopke et al.l2007l) o r dela yed detonations dGamezo et al 



have focused on turbule nt deflagrations dSchmidt & Niemevei The opposite problem occurs in the case of rapidly rotat 



Ropke & Niem everll2007 ) of Chandrasekhar-mass WDs 



ing WDs, which burn as turbulent deflagrations dPfannes et al 



(this class includes the gravita ti onally confined detonation 



model proposed by iPlewa et al.l ([2004)). The delayed deto 



nation model yields an explosion consistent with kinemat- 
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d2009h . tereafter referred to as Paper I). Here, too much stellar 
material remains unburned as a consequence of the lower den- 
sity at large radii. Combining the fast propagation of a deto- 
nation front with the shallower stratification of the stellar fuel, 
rapidly rotating initial models have already been proposed as 
a possible means of repairing the prompt detonation mecha- 
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nism by ISteinmetz et al' ( 1992 ). However, in their models the 
bulk of WD matter was nevertheless burnt to form iron-group 
elements (IGEs) and, accordingly, even rapid rotation did not 
change the situation. Whenever a significant amount of IMEs 
was produced by rotation in their study, the amount of IGEs 
was simultaneously too high. 

In this work, we revisit the scenario of pr omptly detonat- 



ing ca rbon-oxygen (CO) WDs. In contrast to ISteinmetz et al 
(1992), we investigate a variety of rotation laws includ 



ing d ifferentially rotating models inspired by the results of 
lYoon & Langei ( 2005 ). Some of the models in our sample sig- 
nificantly exceed the Chandrasekhar mass , which is indeed sug - 
gested by some observations, for instance. lHowell et al.1 (|2006). 
For this reason, we do not attempt to interpret our results as 
those for normal SNe la but instead look for possible connec- 
tions with observed peculiar supernovae. 

It is important to emphasize that we only attempt to predict 
the expected range of masses of the produced IMEs and IGEs, 
to determine whether these types of explosions might account 
for rare SN la events. To this end, we first provide two sim- 
ple estimates based on the equilibrium stratification of the WD 
and different choices of the density threshold for burning into 
nuclear statistical equilibrium (NSE). These numbers can be in- 
terpreted as plausible upper and lower bounds on the produced 
IME masses in these events. This interpretation is confirmed 
by the nucleosynthetic yields obtained by post-processing a 
full explosion simulation, which lie in between these bounds. 
Going beyond these estimates would require a parameter study 
of rotation laws and ignition conditions of sufficient resolution 
to capture the detonation front. Given the speculative nature 
of this model, we do not think that this effort is warranted at 
present. 

Based on these simplifying assumptions, we find that a sig- 
nificant amount of IMEs (0.1 to 0.4 Mq) is being produced 
accompanied by a Super-Chandrasekhar mass amount of IGEs 
( 1 .5 to 1.8 Mq) and a marginal amount of unburnt stellar mate- 
rial. The ejecta expand at higher radial velocities because of the 
greater amount of nuclear energy released as a consequence of 
the detonation. Unburnt stellar material is leftover only in stel- 
lar regions close to the stellar surface at the equatorial plane, 
whereas iron group elements are predominant within the star 
and at the stellar poles as a consequence of the density stratifi- 
cation that is shaped by rapid rotation. For the same reason, a 
bulge of IMEs is generated in the outer regions of the equatorial 
plane. 

The discussion will proceed as follows. In Sect. 2, we in- 
troduce our numerical method. The procedure that allows us 
to estimate the composition of burning species is discussed in 
Sect. 3. The process of a typical prompt detonation initiated 
within a rapidly rotating WD is described in Sect. 4. Section 5 
contains an interpretation with respect to spectral features. The 
detailed stellar composition is investigated in Sect. 6. Section 7 
concludes the paper. 



2. Method 

To study the outcome of detonations of rapidly rotating WDs, 
we constructed rotators with different rotation laws that can po- 
tentially give rise to SN la progenitors. 

The WDs were assumed to consist of carbon and oxygen 
in equal proportions with central densities of 2.0 x 10 9 g/cm 3 . 
Depending on the detailed rotation law, their masses ranged 
from 1.4 to 2.1 Mq (cf. Fig. [2J. The equation of state consid- 
ered the thermodynamic al properties of a photon, electron, and 
baryonic gas, as well as electron-positron-pair-creation. The 
temperature profile was adjusted to describe the evolut ion of 
the accreting WD (cf. the study of lYoon & Langerl l 20051) ). Our 
simulations were carried out in three spatial dimensions. We 
used a moving equidistant cartesian grid with 100 3 cells (ini- 
tial grid spacing of 1.0 x 10 7 cm) for the computation of a full 
star explosion. The detonation was ignited centrally for most of 
the calculations. Gravitational effects were taken into account 
by using a multipole solver (see Paper I for more details of the 
rotation laws and the numerical implementation). 

Since the detonation front incinerates the star before 
changes in the stellar structure can occur, an estimate of the 
burning products of detonating WDs can be obtained by sim- 
ply looking at the material densities present in the hydrostatic 
case. 

The thermal and chemical structu re of planar thermonu- 
clear detonations is well understood dlmshennik & Khokhlov 
19841) and yields a transition density for burning into NSE of 
approximately 10 7 g/cm 3 . However, in the models that we con- 
sider, the detonation front is subject to strong s hear, which, in 



combination with the cellular front instability (iGamezo et al 



1999 r.lTimmes et alJ2000h can be expected to give rise to strong 
local variations in the temperature and density profile. Most of 
these effects have a tendency to lower the burning tempera- 
ture at fixed density and thus postpone the burning to NSE to 
a higher threshold. A full exploration of these phenomena is 
beyond the scope of our present study. Instead, we argue that 
the higher transition density of planar deflagrations (roughly 5 
xlO 7 g/cm 3 ) yields a plausible upper cutoff, the lower cutoff 
being given by planar detonations. 

Hence, to estimate the range of masses of IMEs produced 
by the detonation, we considered two different burning thresh- 
olds for NSE that can be considered as limiting cases. As an 
upper bound to the transition density, we applied the same 
density thresholds as commonly used for deflagrations. There, 
within regions of high density (p > 5.248 x 10 7 g/cm 3 ), 
IGEs are assumed to be produced. Medium density regions 
(5.248 x 10 7 g/cm 3 > p > 1.047 x 10 7 g/cm 3 ) and low den- 
sity regions (p < 1.047 x 10 7 g/cm 3 ) yield IMEs and unburnt 
fuel ("C+O"), respectively. We refer to those high NSE burning 
density threshold values hereafter as the HBT limiting case. 

The transition density for planar detonations was chosen to 
be a lower bound (low NSE burning density threshold, in short 
LBT). IGEs are assumed to be produced for p > 10 7 g/cm 3 , 
IMEs are obtained for lower densities. 

Finally, to take into account the explosion dynamics, hy- 
drodynamical simulations are necessary. Thus, in addition to 
the analysis of the density stratification, we used the hydrody- 
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Fig. 1. Fractions of the total mass in the high density regime 
(p > 5.248 x 10 7 g/cm 3 , "IGEs"), the medium density regime 
(5.248 x 10 7 g/cm 3 > p > 1.047 x 10 7 g/cm 3 , "IMEs"), and 
the low density regime (p < 1.047 x 10 7 g/cm 3 , "C+O") for 
different rotation laws versus the ratio of equatorial to polar 
radius (denotation according to the HBT limiting case) 

namics code PROMETHEUS (cf. Paper I) to simulate prompt 
detonation s of rotating WDs. By employing the method de- 
scribed by iTravaglio et al. ( 2004) . the data from these simu- 
lations were post-processed to derive a more accurate result for 
the explosive nucleosynthesis. We emphasize, however, that the 
resolution and range of explored parameters are insufficient to 
provide a robust investigation of the nuclear yields of this sce- 
nario. 

Burning i s realized by me a ns of the level set method, in- 

and applied to SN la 
for pure deflagrations. 



traduced by Osher & Sethian 



calculations by Reinecke et al 




Golombek & Niemever (2005) applied this method to super- 
sonic burning fronts in the context of delayed detonations. In 
both cases, the flame front is treated as a surface that separates 
fuel from ashes with a sharp discontinuity. For detonations, the 
speed of the flame depends on the mass density and i s calcu- 
lated by a linear interpolation of the values in Fig. 2 of lSharpe 
dl999h . For the hydrodynamical simulations we used the HBT 
limiting case, taking into account that temperature is underes- 
timated once densities decrease below 5.248 x 10 7 g/cm 3 . 

3. Estimate of the burning products 

Owing to detonations being able to propagate through an es- 
sentially undisturbed equilibrium configuration, we can esti- 
mate the burning species on the basis of the hydrostatic density 
stratification of the progenitor star. 

While the detonation front crosses dense material, 56 Ni is 
mainly produced (for densities p > 3 x 10 9 g/cm 3 , although 
adverse electron capturing becomes important). A sufficient 
amount of IMEs can only be obtained if there is enough mate- 
rial at densities lower than about 5xl0 7 g/cm 3 or 1 xlO 7 g/cm 3 
(HBT and LBT limiting cases, respectively) in the progeni- 
tor star. Figure Q] illustrates the mass proportions in the corre- 
sponding density regimes for the rigid rotator, a rotator with ap- 
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Table 1. Quantities for the "AWD3 detonation" explosion 
model from the initial model in the HBT and LBT limiting 
cases. 



proximately constant angular momentum (j con si), and a rotator 
obeying the rotation law derived for an accreting white dwarf 
(AWD; cf. Paper I and the study of lYoon & Langerld2004 ). For 
every kind of rotator, the ratio of the equatorial radius to the po- 
lar radius steadily increases along the x-coordinate and the WD 
is deformed into a doughnut-like shape with increasingly rapid 
rotation. 

From Fig.Q] it is clear that for the AWD3 series, fewer IGEs 
and, at the same time, more IMEs are generated in comparison 
to the other models. The AWD3 rotator with the fraction of the 



radii r equator /r po i = 2.4 contains 1.41 Mq of IGEs and 0.45 Mq 
of IMEs, or 1.86 M Q of IGEs and 0.21 M Q of IMEs (HBT 
and LBT limitin g cases, respectively). In contrast, the j CO nst se- 
quence used by ISteinmetz et al] (119921) produces more IGEs 
and fewer IMEs. 



4. Evolution of a prompt detonation 

Figure [2] shows the density contours of the centrally ig- 
nited "AWD3 detonation" explosion model. The burning front 
(white, thick line) moves outwards (without stopping after 
burning has ceased but continues to propagate outside the star 
for numerical reasons). The large amount of IMEs resulting 
from the prompt detonation of the AWD3 rotator stems from 
the equatorial bulge of low density material that is typical of ac- 
creting WD rotation (the density thresholds for the production 
of IMEs and IGEs, respectively, are indicated by the contour 
lines inside the burning front for the HBT limiting case). 

Table Q] summarizes the values derived directly from the 
rotating initial model for both HBT and LBT limiting cases 
(left and middle column respectively). The ratio of rotational 
to gravitational energy of the progenitor star is denoted by /3. 
The energies are obtained from the hydrodynamical simula- 
tion. Moreover, the results from post-processing of the simu- 
lation data are listed (cf. Sect. [5). Remarkably, the results of 
the post-processed simulation are marked on either side by the 
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Fig. 2. Density contour lines for the prompt detonation of the AWD3 rotator at different instants of time (the white, thick line 
represents the detonation front). Cross-sections along the rotation axis of a full star simulation are shown. 



limiting case estimates of the equilibrium model, confirming 
the consistency of our assumptions. 

The mass of the tracer particl es used in post-processing the 



simulation data dReinec ke 2001) adds up to 2.12 Mq, which 



exceeds the stellar mass of the AWD3 rotator by more than 2 %. 
This is a consequence of the small number of tracer particles, 
namely 15 3 , employed by the post-processing, which is com- 
putationally demanding because of the high detonation temper- 
atures. 

Compared to the subsonic deflagrations, the entire explo- 
sion proceeds faster for detonations because of the supersonic 
propagation speed. Figure [3] shows the evolution of the energy 
contributions of the turbulent deflagration of a similar rotator 
and the "AWD3 detonation" explosion model. The high tem- 
peratures caused by the detonation lead to a noticeable peak in 
the internal energy at t ~ 0.1 s. The nuclear energy released 
by the detonation exceeds E nuc for the deflagration scenario by 
the factor of ~ 2.5. This energy release is sufficient for an early 
unbinding of the star (which is only marginally achieved in the 
deflagration case for the specific deflagration scenario) and the 



ejection of material with a kinetic energy that is increased by a 
factor of ~ 3. 

As previously mentioned, the carbon-to-oxygen ratio was 
fixed to X( 12 C) = 0.5 for the detonation study. However, sev- 
eral publi cations indicate that in the WD's interior it is less than 
0.5 dCouch & ArnetJl975l;IUmeda et al.ll999tlHoflich & Stein 



l2002t iLesaffre et al.ll2006l) . which would produce a less ener- 
getic explosion. 



5. Expected spectral features 

The most striking property with respect to spectral features is 
the extreme distribution of species in the homologous expan- 
sion phase. Figure [4] shows cross-sections along the rotation 
axis for the total and fractional mass densities at t — 2.3 s. 
The spatial extent of the ejecta at that time corresponds ap- 
proximately to the expansion reached for pure deflagrations af- 
ter t — 5 s (cf. Paper I). Whereas IMEs are present only in 
a torus within the equatorial plane, IGEs emerge close to the 
poles and are therefore already visible at an early period. No 
fuel is left at the centre of the star. The spatial distribution of 
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Fig. 4. Total (upper left) as well as fractional mass densities (upper right: IGEs; lower left: unburnt fuel; lower right: IMEs) of 
the AWD3 rotator after t = 2.3 s (HBT limiting case). 



the nuclear species should produce a strong dependence on the 
line of sight. 

Figure [5] shows the resulting radial velocities at t = 5 s for 
the different species in comparison to the deflagration model. 
No fuel appears at low velocities, which is a necessary fea- 
ture of an SN la model. IGEs are encountered over a broad 
range of velocities but are most prominent at lower radial ve- 
locities (v r ~ 5 x 10 3 km/s). IMEs are found at higher ve- 
locities, and are most prominent at v r ~ 10 x 10 3 km/s. In 
addition, some of the ejecta exceed velocities of up to v r ~ 
20 x 10 3 km/s. The refore, high velocity features (HVFs; cf. 
Mazz ali et al. [ (Eoolh could arise from prompt detonations of 
rapidly rotating WDs. The composition structure and velocity 
satisfy co nstraints that deriv ed for SN 199 IT using a synthetic 
spectrum ( Fisher et al . 1999) . Additional benefits of the prompt 
detonation of the AWD3 rotator are given by wor k on the three- 
dime nsional spectral synthesis for early spectra (ITanaka et alj 
2006). It was found that a thick torus — similar to the one em- 
bodied by the IMEs — can naturally explain the diversity in 
observed strength of the HVFs by the covering of parts of the 
star. 



6. Post-processing of the prompt detonation 

By considering a detailed nuclear reaction network, we gained 
insight into the composition resulting from the "AWD3 detona 



tion" s cenario by means of post-processing (cf. lTravaglio et al 



(2004) 'or details of the method). In particular, we compared 
the composition obtained from the post-processing to the com- 
position resulting from the density-dependent burning of the 
hydrodynamical simulation. Even though differences in the de- 
rived compositions are discernable (see Table [T] in which all 
species exceeding the amount of 5 x 10~ 3 Mq are listed), a 
significant amount of IMEs resulting from the prompt detona- 
tion of an accreting WD rotator is produced. 

The detailed composition calculated by means of post- 
processing is listed in Table|2] The total amount 0.29 Mq gen- 
erated within the AWD3 rotator compares fairly well to the SN 
1 99 1 T-like events, which exhibit more than a solar mass of 56 Ni 
and a small amount of 0.2 to 0.3 Mq of IMEs (P. Mazzali, pri- 
vate communication). Table|2]lists all the isotopes of IGEs and 
IMEs that exceed five-tenths of a percent of a solar mass. The 
amount of 1 .48 Mq of 56 Ni is extraordinarily high and should 
lead to very bright SNe la. As a result of high temperatures, 
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Table 2. Composition of various species for explosion scenario 
of "AWD3 detonation" from the post-processing calculation. 



fewer neutron-rich isotopes among IGEs are produced. For in- 
stance, the ratio of 56 Ni / 58 Ni — 11. Whereas the radioactive 
decay of 56 Ni supplies the energy for the light emitted by the 
SN la ejecta, the presence of other iron-group nuclei enhances 
the op acity of the ejecta an d causes a broadening of the light 
curve (Reinecke et al ]l2002h . 




2.0-10' 



(b) "AWD3 detonation" 

Fig. 5. Probability density functions in radial velocity space 
for the deflagration scenario (a) and the prompt detonation 
(HBT limiting case) (b) at t — 5 s. The velocity increments 
are different for the deflagration and detonation simulation 
(dv = 5 x 10 6 cm/s and dv = 4 x 10 7 cm/s, respectively) 
as a result of the different resolution for both hydrodynamics 
simulations. 



7. Detonating rotators: only an exotic scenario? 

We have demonstrated that the prompt detonation of a rapidly 
rotating WD potentially meets the criteria of superluminous 
SNe la regarding spectral and bolometric features. In par- 
ticular, prompt detonations are capable of producing a suf- 
ficient amount of IMEs, if the progenitor star is a super- 
Ch andrasekhar-mass WD with differential rotation as proposed 
by lYoon & Langerl d2005h . Therefore, this explosion scenario 
may provide an explanation of the rather small number of 
events such as SN 199 IT or SN 2003fg. 

In our investigation, we assumed that prompt detonations 
can occur in rapidly rotating WDs. For WDs that do not ro- 
tate appreciably, on the other hand, pure deflagration or de- 
layed detonation models, where the DDT occurs at a relatively 
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late stage of the explosion, are most successful. The prompt 
detonation may be seen as a limiting case of a defiagration-to- 
detonation transition, where the transition time approaches zero 
as the rotation frequency of the WD becomes nearly critical. 
But this appears to be implausible, if the mechanis m causing 



DDT s is related to the int ermittency of turbulence dPan et al 



2008; Schmid t et al.l 120091) . Alternatively, prompt detonations 
might occur beyond a certain rotation frequency, whereas deto- 
nations are delayed or entirely absent for WDs rotating at lower 
frequencies. However, given the rudimentary theoretical un- 
derstanding of the onset of explosive thermonuclear burning 
in WDs (both non-rotating and rotating), no mechanism that 
selects the initial burning mode depending on the rotation fre- 
quency is known at present. 

Acknowledgements. We thank Fritz Ropke for post-processing the 
"AWD3 detonation" explosion model, Paolo Mazzali for discussing 
the results, and the referee for helpful remarks about the NSE transi- 
tion densities. 
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